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We find the most general spherically symmetric non singular black hole solution in a 
special class of teleparallel theory of gravitation. If r is large enough, the general solution 
coincides with the Schwarzschild solution. Whereas, if r is small, the general solution be- 
^ | haves in a manner similar to that of de Sitter solution. Otherwise it describes a spherically 
symmetric black hole singularity free everywhere. Moreover, the energy associated with the 
^ 'general solution is calculated using the superpotential given by M0ller 1978. 
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1. Introduction 



M0ller modified general relativity by constructing a new field theory in teleparallel space 
[jy . The aim of this theory was to overcome the problem of energy-momentum complex that 
appears in the Riemannian space 0. The field equations in this new theory were derived 
from a Lagrangian which is not invariant under local tetrad rotation. Sdez || generalized 
M0ller theory into a scalar tetrad theory of gravitation. Meyer jj| showed that M0ller theory 
is a special case of Poincare gauge theory [|, 0. 

Mikhail et al. [0 obtained two spherically symmetric solutions in M0ller's theory when 
the stress-energy momentum tensor vanishing. They calculated the energy associated with 
those solutions and found that in one solution the energy does not coincide with the gravita- 
tional mass. Shirafuji et al. f| extended the calculation to all the stationary asymptotic flat 
space solutions with spherical symmetry, dividing them into two classes, the first in which 
the components (e a0 ) and (e 0Q ) of the parallel vector fields (e^) tend to zero faster than 
1 / y/r for large r and the second, in which these components go to zero as 1/ ^fr. It was 
found that the equality of the gravitational and inertial masses holds only in the first class. 
Mikhail et al. M also obtained a spherically symmetric solution of M0ller's theory starting 
from a tetrad which contains three unknown functions following Mazumder and Ray |Uj. 
The solution contains one arbitrary function of the radial coordinates r, and all the previous 
solutions can be obtained from it. 

Dymnikova derived a static spherically symmetric nonsingular black hole solution in 



orthodox general relativity assuming a specific form of the stress-energy momentum tensor. 
This solution practically coincides with the Schwarzschild solution for large r, for small 
r it behaves like the de Sitter solution and describes a spherically symmetric black hole 



singularity free everywhere It has been proved that it is possible to treat this specific 
form of the stress-energy momentum tensor as corresponding to an r-dependent cosmological 
term A Mi/ , varying from A^ v = Ag^ u as r — > to A^ v = as r — ► oo [f[2| . More recently 



131 , the spherically symmetric nonsingular black hole has been used to prove that a baby 
universe inside a A black hole can be obtained in the case of an eternal black hole. Also 
it has been shown that the probability of a quantum birth of a baby universe can not be 
neglected due to the existence of an infinite number of A white hole structures. Assuming 



the specific form of the vacuum stress-energy momentum tensor using a spherically 



symmetric tetrad contains three unknown functions of radial coordinate r only, which was 



constructed by Robertson WM , the author |T5| obtained a general solution of this tetrad when 



the stress-energy momentum tensor is not vanishing. This solution contains one arbitrary 
function of the radial coordinates r, and all the previous solutions could be derived from it. 

The general form of the tetrad, e^, having spherical symmetry was given by Robertson 



141. In the Cartesian form it can be written a: 



e ° = iA, e a ° = Cx a , e a = iDx a 

e a a = b a a B + Fx a x a + ea^Sx 13 , (1) 

where A, C, D, B, F, and S are functions of t and r = (x*!*) 1 / 2 , and the zero th vector eo M 
has the factor i = ^—1 to preserve Lorentz signature. We consider an asymptotically flat 

*In this paper Latin indices ...) represent the vector number, and Greek indices (/i, v, ...) represent 
the vector components. All indices run from to 3. The spatial part of Latin indices are denoted by (a, b, ...), 
while that of Greek indices by (a, (3, ...). In the present convention, Latin indices are never raised. 



space-time in this paper, and impose the boundary condition that for r — > oo the tetrad (1) 
approaches the tetrad of Minkowski space-time, (e^) = dia.g(i,5 a a ). 

The aim of the present work is to find a general solution with spherical symmetry in 
M0ller's tetrad theory of gravitation when the vacuum stress-energy momentum tensor is 



not vanishing and has the same form given by Dymnikova |]TlJ . Moreover, the energy of that 
solution is calculated. In section 2 we briefly review M0ller's tetrad theory of gravitation. In 
section 3, we study the general solution without the S-term (see 1), where the remaining un- 
known functions are allowed to depend on t and r and the general solution with an arbitrary 
function of t and r is obtained. Then we compare this general solution with that obtained 



before fllpl , it is found that this general tetrad is just the t- independent case of our general 
tetrad without the S'-term. Also we first study the general, spherically symmetric solution 
with a non-vanishing S'-term and a solution with one parameter is obtained. In section 4 
we calculate the energy content of those general solutions using the superpotential given by 
M0ller in 1978. The final section is devoted to discussion and conclusion. 

Computer algebra system Maple 6 is used in some calculations. 



2. M0ller's tetrad theory of gravitation 



In a space-time with absolute parallelism the parallel vector fields define the nonsym- 
metric connection 

F A d = e V (2) 

where = d^e^. The curvature tensor defined by T x ^ is identically vanishing, however. 

M0ller's constructed a gravitational theory based on this space-time. In this theory the 
field variables are the 16 tetrad components e^, from which the metric tensor is derived by 

ST A = c* V. (3) 

We assume an imaginary values for the vector in order to have a Lorentz signature. We 
note that, associated with any tetrad field there is a metric field defined uniquely by (3), 
while a given metric g^ v does not determine the tetrad field completely; for any local Lorentz 
transformation of the tetrads leads to a new set of tetrads which also satisfy (3). The 
Lagrangian L is an invariant constructed from 7 W and g^, where 7 Mi/ p is the contorsion 
tensor given by 

Ipvp ('i p&iu; pi (4) 

where the semicolon denotes covariant differentiation with respect to Christoffel symbols. 
The most general Lagrangian density invariant under the parity operation is given by the 
form 



where 



9 d = det(g^), (6) 



and $ p is the basic vector field defined by 



^ d = TV (7) 

Here a±, 012, and are constants determined by M0ller such that the theory coincides with 
general relativity in the weak fields: 

ai = --, a 2 = -, a 3 = -(l-2A), (8) 

Kj Kj K 

where k is the Einstein constant and A is a free dimensionless parameterQ. The same choice 
of the parameters was also obtained by Hayashi and Nakano . 

M0ller applied the action principle to the Lagrangian density (5) and obtained the field 
equation in the form 

G^v + Hp,v = -kT^, (9) 
Ffiv = 0, (10) 
where the Einstein tensor G^ v is defined by 

G(w — — -^9^vR- (11) 



Here H„ u and F„ v are given by 



u def. , 



(12) 



and 

F, u d = A - - $ p ( 7 V - l P v») + T„/J , (13) 
and they are symmetric and skew symmetric tensors, respectively. 

M0ller assumed that the energy-momentum tensor of matter fields is symmetric. In 
the Hayashi- Nakano theory, however, the energy- momentum tensor of spin- 1/2 fundamental 
particles has non-vanishing antisymmetric part arising from the effects due to intrinsic spin, 
and the right-hand side of (10) does not vanish when we take into account the possible effects 
of intrinsic spin. 

It can be shown |17| that the tensors, and F^, consist of only those terms which are 
linear or quadratic in the axial- vector part of the torsion tensor, a M , defined by 



a, d = \^ P «l vp \ (ID 



where e^ upa is defined by 

tfivpir ~ \J 9^ pupa i (15) 

where 5 M ^ po - being completely antisymmetric and normalized as #0123 = — 1- Therefore, both 
and F pv vanish if the a M is vanishing. In other words, when the a M is found to vanish 
from the antisymmetric part of the field equations, (10), the symmetric part (9) coincides 
with the Einstein equation. 



"Throughout this paper we use the relativistic units, c = G = 1 and k = 8tt. 



3. Spherically symmetric solutions 



In this section we find the most general, spherically symmetric non singular black hole 
solution of the form (1) in Miller's theory. The axial-vector part of the torsion tensor, a M , 
is vanishing, and the skew part of the field equation is satisfied identically as is explained 
above. We discuss two cases separately: One with S = and the other with S ^ 0. 



(i) The case without the S'-term. In this case the axial-vector part of the torsion tensor 



a M is identically vanishing. Thus, when this tetrad is applied to the field equations, the skew 
part (10) is automatically satisfied and the solution of the symmetric part is the nonsingular 
black hole solution given before by Dymnikova [IT|. Therefore, the solution of the form (1) 
with S = can be obtained from the diagonal tetrad of the nonsingular black hole metric 
by a local Lorentz transformation which keeps spherical symmetry, 



(A 



kl) 



( L iH sin 9 cos iH sin 9 sin 

— iH sin9 cos <p 1 + (L — 1) sin^ 2 cos0 2 (L — 1) sin# 2 sin cos < 

— iH sm. 9 smip {L — 1) sin # 2 sin cos 1 + (L — 1) sin 9 2 sin 

\ — iH cos 9 (L — 1) sin # cos # cos (L — 1) sin # cos # sin </ 

where H is an arbitrary function of t and R and 



iH cos 9 

(L — 1) sin 9 cos 9 cos < 

(L — 1) sin 9 cos 9 sin < 

1 + (L - l)cos# 2 
(16) 



L = VJP + l. 

e/ = W 0)AI 



Namely, we see that 

(17) 

is the most general, spherically symmetric solution without the S'-term. Here e/ (0)/i is the 
diagonal tetrad in the spherical polar coordinates given by [|18[ 
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and R is defined as R = r/B (T^j. The explicit form of the is then given by 
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where T\ is defined by 
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(19) 
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(20) 



If we apply the tetrad (19) to the symmetric part of the field equation (9), the right hand 
side takes the form 



T ° = Tl 1 = e e 



T 2 2 



T,3 



e e 



-R 3 /ri 



3R* 
2r! 3 



(21) 



The metric associated with the tetrad (19) is given by 

ds 2 = -£dt 2 + 



2 < ^ ' dR2 +R 2 dn 2 , 



where 



i = X< 



and 



dtt 2 = d9 2 + sin 2 



:22i 



(23) 



which is the spherically symmetric nonsingular black hole solution [11 



(ii) The case with non-vanishing S-term. We start with the tetrad of (1) with the six 



unknown functions of t and r. In order to study the condition that the vanishes it is 
convenient to start from the general expression for the covariant components of the tetrad, 



eoo 



iA, e a o = Cx a , eoa = iDx c 
5 aa B + Fx a x a + e aaf3 Sx^ 



(24) 



where the six unknown functions, A, C, D, B, F, S, are connected with the six unknown 
functions of (1). We can assume without loss of generality that the two functions, D and F, 
are vanishing by making use of the freedom to redefine t and r [17]. We then transform the 
tetrad (24) to the spherical polar coordinates (r, 9, 0, t) to take the form 



/ iA 
rC sin 9 cos < 
rC sin 9 sin ( 
V tC cos 9 





B sin 9 cos < 
B sin 9 sin < 
Bcos9 





rB cos 9 cos + r 2 S sin i 
tB cos 9 sin — r 2 S cos < 
—rB sin 9 





—rB sin 9 sin + r 2 S sin 9 cos 9 cos < 

rB sin 9 cos + r 2 S sin 9 cos 9 sin 

-r 2 5 sin 2 9 
(25) 



Then the condition that the axial vector part a M vanishes is to be || . 



n fTTi v \ 3BS + r< - BS ' - - S ' S '- f = °' 

° = V<-9»°'=| 2CS + (SB-SB), „=! (26) 

with S" = dS/dr and S = dS/dt. This condition can be solved to give 

(7 = 0, 5 = ^B, (27) 

where rj is a constant with dimension of (length) 2 . 

The symmetric part of the field equations now coincides with the Einstein equation. The 
metric tensor formed of the tetrad (25) with (27) is not of the non singular black hole solution 
form. Taking the new radial coordinate 



R = rBJl + t (28) 



then the metric tensor takes the well-known non singular black hole solution form [ IT ]. 

Applying the transformation (28) to the tetrad (25) with (27), we finally obtain the 
general, spherically symmetric non singular black hole solution with non-vanishing S'-term: 
The non-vanishing, covariant components of the tetrad are given by 

e 00 = iX 

sin 9 cos 6 



en 



X 

R 

ei2 = — {a r] f 2 (R) sin + 4 f 4 (R) cos 9 cos 4> — rj 2 cos 9 cos 
R 

ei3 = — (4 77 f 2 (R) cos 9 cos — 4 / 4 (i?) sin <p + rf sin 0J sin 9 
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sin 9 sin 1 
X~~ 



622 = —y- ^4 r] f 2 (R) cos0 — A f i (R) cos # sin <\> + r/ 2 cos sin 
e 23 = 77 (4 77 f 2 {R) cos 6* sin + 4 f A (R) cos <p — f] 2 cos sin # 



y 

cos 



e 3 i 



e 32 = ^(4/ 4 (i?)-r ? 2 ) sin# 

es 3 = ^vf 2 (R) sin 2 #, (29) 



where 



X = Jl - ^ (l - e- rf M 8 ), y = 4/ 4 (i?) + r, 2 



and /(-R) is given by 

/( J R)=e^ dR /( RX ), (30) 



The tetrad (29) when applies to the symmetric part of field equation (9), the right hand 
side takes the form 



T ° = Ti 1 = e e 



-R 3 / ri 



T 2 2 = T 3 3 = e e~ R3 ^ I 1 



The metric associated with the tetrad (29) is given by 



3R 3 
2ri 3 



(31) 



jd2 

1 < 1 ' + R 2 dn 2 , 



ds z = -^dt 2 + 



(32) 



where 



(33) 



which is the spherically symmetric nonsingular black hole solution |11 



It is clear that if rj and H (R, t) are equal to zero the two classes of solutions given by 
(19) and (29) coincide with each other, and reduce to the solution given before | 18| . Fur- 
thermore, if the exponential term is equal to zero then the two solutions coincides with each 
other and give the Schwarzschild solution. 

The solutions (19) and (29) are the exact solutions of the M0ller's field equations. They 
practically coincide with the Schwarzschild solution for R » r\ or R >> r\ and behave 
like the de Sitter solution, for R « r\ or R « r\. 

As is clear from (21), the spherically symmetric stress-energy momentum tensor is really 
anisotropic. The difference between the principle pressures 



-Pk, 



(34) 



correspond to the well known anisotropic character of evolution of the space-time inside 
a black hole undergoing a spherically symmetric gravitational collapse [ 19| . For R << r\ 



isotropization occurs and the stress-energy momentum tensor takes the isotropic form 



T a /3 — eg a/3. 



(35) 



When R — ► the energy density tends to eo- For R>> r± all the components of the stress- 
energy momentum tensor tend to zero very rapidly. The same properties satisfied for the 
stress-energy momentum tensor (31). 

Now let us compare the solution (19) with that given before ||15|| . Started from a spheri- 
cally symmetric tetrad with three unknown functions of the radial coordinate r only, which 
is given in the spherical polar coordinates by 
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(36) 



Applying the tetrad (36) to the field equations, (9) and (10), we obtained a general solution 
of the form 



A(R) 
V(R) 
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1 [2m 
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(37) 



where A(R), B(R) and D>(R) are connected to A, B and D through the transformation 

R = r/B, B' = — -5—, and B is an arbitrary function of R. The tetrad (36) can further 
dR 



be expressed in the form [15 
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/ 



Thus, it is easy to verify that this tetrad can be obtained from (19) by choosing the function 
H as 



H 



R 2 B' 2 -2RB' + ^(l-^- k3 ^ 

r_A 



1/2 



2m 
R 



(39) 



4. The Energy Associated with each Solution 



The superpotential of the M0ller's tetrad theory of gravitation is given by Mikhail et al. 
as 

-g) 1/2 



2k 



-p xp J vX [® p g ax g„T - *9r»l xpa - (i - 2X)g Tlll 



where P XP J vX is 



p TvX d = S T a vX 4- S T a uX - S r a uX 



(40) 
(41) 



with g pa u being a tensor defined by 



(42) 



The energy is expressed by the surface integral 



20 



E = lim / U 0a n a dS, (43) 

r ^°° J r=constant 

where n a is the unit 3-vector normal to the surface element dS. 

Now we are in a position to calculate the energy associated with the two solutions (19) 
and (29) using the superpotential (40). As is clear from (43), the only component which 
contributes to the energy is U§ a . Thus substituting from the solution (19) into (40) we 
obtain the following non-vanishing value 

= ^C(L-X). (44) 
Substituting from (44) into (43) we get 

E(R) =Xh(L-X). (45) 

As is clear from (45) that the energy depends on the arbitrary function H. If this arbitrary 
function takes the value (39) then 

E(R) = 2m (l - e -^ 3/ri3 ) - R 2 B', (46) 



which is the same as that obtained before [T^ . 



Now let us turn our attention to the solution (29). Calculating the necessary components 
of the superpotential, we get 



u " — r y - x )- (47) 

Substituting from (47) into (43) we get 

E(R) = XR^%-" 2 >- X y (48) 

where f(R) is given by (30). As is clear from (48) that if the constant rj — 0, then 

E(R) =m(l-e- R3/rii ), (49) 



which is the same as that given before fll8[ , Yang |2T| and Radinschi ||22|| . This is a very 
satisfactory result since the two tetrads (29), (36) coincide when r] = and D = 0, and 
accordingly their energies must be identical. We can also see from (48) that if R — *■ 0, then 
E(R) —>■ and form (49) if R — » oo then E(R) — > m. Thus the total energy will coincide 
with that of the Schwarzschild solution. 



5. Main results and Discussion 



We have studied the general spherically symmetric nonsingular black hole solutions in 
teleparallel theory of gravitation. Tetrads with spherical symmetry are classified into two 
groups. 

i) The groups in which the S-term is vanishing, and therefore, the axial vector part is 
identically vanishing. Accordingly we obtain a family of solutions with arbitrary function of 
t and R, i.e., H(t, R). A relation with the solution given before [ T5| had been established by 



assuming a specific form of the arbitrary function H(t, R) as given by (39). 

ii) The groups which possess the S-term, namely, the term S(t, r)e aa px^ . When the S-term 
is non-vanishing, we obtain the condition (27) which is the condition for vanishing the axial 
vector part, i.e. a M = 0. In this case we obtain a family of solutions with a constant 
parameter, i.e., rj. For comparison if this constant is set equal to zero then the tetrad (29) 
will reduce to the diagonal form of (36). 

Pereira et al. []23| studied the teleparallel version of the exact solutions ( Schwarzschild 
and Kerr solutions) of general relativity. As is well known that the torsion tensor can be 



decomposed into three parts tracless part, vector part and the axial vector part [fL7l| . Using 



the weak-field limit they |23| show that the vector and the tracless parts of the Schwarzschild 



torsion combine to yield the Newtonian force. By considering the slow rotation and weak- 
field approximation they show that for the Kerr solution the axial vector part is nothing but 
the gravitomagnetic components of the gravitational field and is therefore the responsible 
for the Lense Thirring effect. 

We also calculated the energy associated with the two solutions (19) and (29). Concerning 
the energy associated with the solution without the S-term, it is found that the energy 
depends on the arbitrary function H(t,R). If the function H(t,R) takes the value (39), 



then the energy coincides with that obtained before |15] and if H (t, R) = 0, then the energy 
coincides with that obtained by Yang [^T[] and Radinschi 

As for the energy associated with the solution when the S-term is non-vanishing we 
obtain the expression (48). It is clear from this expression that the energy depends on the 
parameter 77 and if this parameter equals zero then the tetrad (29) will coincide with the 
diagonal form of the tetrad (36) and the associated energy will be the same as that obtained 
before 0. It is clear from (48), that if R -> 0, then E(R) -> and from (49) if R -> 00 
then E(R) — > m. Thus the total energy will coincide with that of the Schwarzschild solution, 
since the metric associated with the tetrad (29) in this case gives the Schwarzschild metric. 

It is readily seen from equations (46) and (49) that the factors in front of the gravitational 
mass m are different. This is because that when H(t, R) = 0, the components (e a o) and (eo a ) 

of the parallel vector fields (19) tend to zero as l/\f~R\ and in that case the gravitational 
mass will not coincide with the inertial mass even in the simple case of vacuum, i.e., eo = 
[0. ||, 3. In contrast to the tetrad (29) when r\ = 0, the components (e a o) and (eo Q ) tend 
to zero faster than 1 / y/~R and this case is satisfactory since the gravitational mass equals to 
the inertial mass f7|, §. 

In the case of the ADM mass we obtain the same result as we can see from (49) 

M ADM = E(R) R ^ oc = m. 

But for equation (46) the result will depend on the form of the arbitrary function B. 

We calculate the energy using expression (43) mainly given by M0ller [|IJ. Blagojevie' and 
Vasilie' |25[ presented an investigation of the connection between the asymptotic Poincare' 



symmetry of space-time and the related conservation laws of energy-momentum and angular 
momentum in teleparallel theory of gravity They |25j derived the generators of the global 
Poincare' symmetry in the asymptotic region from the related gauge generators. They |25 



concluded that the Poincare' generators have to be improved by adding certain surface terms, 
which represent the values of energy-momentum and angular momentum of the gravitating 
system. These results for energy-momentum and angular momentum are valid for general 
teleparallel theory of gravity. 

Chang et al. p6| argued that every energy-momentum complex is associated with a le- 
gitimate Hamiltonian boundary term, and, because of this the energy-momentum complexes 
are quasi-local and acceptable. Each is the energy-momentum density for some physical 
situation. This Hamiltonian approach to quasi-local energy-momentum rehabilitates of the 
energy-momentum complexes. 

A summary of the main results is given in the table below. The general solution of spher- 
ically symmetric nonsingular black hole is classified into two groups according to whether or 
not the space-space components e a a have the term Se aa) 3X^ (referred to as the S-term for 
short). The general solution without the S'-term has an arbitrary function of t and r. The 
general solution with the .S'-term has a constant parameter. 
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